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14.  ABSTRACT 

We  have  adapted  matched  field  processing — a  method  developed  in  underwater  acoustics  to  detect  and  locate  targets —  to  classify  transient  seismic  signals  arising  from  mining 
explosions.  Matched  field  processing,  as  we  apply  it.  is  an  empirical  technique,  using  observations  of  historic  events  to  calibrate  the  amplitude  and  phase  structure  of  wavefields  incident 
upon  an  array  aperture  for  particular  repeating  sources.  The  objective  of  this  project  is  to  determine  how  broadly  applicable  the  method  is  and  to  understand  the  phenomena  that  control 
its  performance.  We  are  inve.stigating  whether  similar  results  can  be  expected  in  different  regions,  with  apertures  of  increasing  scale  and  for  diffu.se  seismicity.  We  verified  .similar 
perfonnance  with  the  closely-.spaced  Zapolyami  mines,  though  discovered  that  it  may  be  necessary  to  divide  event  populations  from  a  single  mine  into  identifiable  subpopulations.  For 
this  purpose,  we  perform  cluster  analysis  using  matched  field  statistics  calculated  on  pairs  of  indiv  idual  events  as  a  distance  metric.  In  our  initial  work,  calibrations  were  derived  from 
ensembles  of  events  ranging  in  number  to  more  than  100.  We  are  considering  the  perfonnance  now  of  matched  field  calibrations  derived  with  many  fewer  events  (even,  as  mentioned, 
individual  events).  Since  these  are  high-variance  estimates,  we  are  te.sting  the  use  of  cross-channel,  multitaper,  spectral  e.stimation  methods  to  reduce  the  variance  of  calibrations  and 
detection  .statistics  derived  from  single-event  ob.servations. 

To  test  the  applicability  of  the  technique  in  a  different  tectonic  region,  we  have  obtained  four  years  of  continuous  data  from  4  Kazakh  arrays  and  are  extracting  large  numbers  of  event 
segments.  Our  initial  results  using  1 32  mining  explosions  recorded  by  the  Makanchi  array  are  similar  to  those  obtained  in  the  European  Arctic.  Matched  field  processing  clearly 
separates  the  explosions  from  three  closely-spaced  mines  located  approximately  400  kilometers  from  the  array,  again  using  waveforms  in  a  band  (6-10  Hz)  normally  considered 
incoherent  for  this  array.  Having  reproduced  ARCES-type  perfonnance  with  another  small  aperture  array,  we  have  two  additional  objectives  for  matched  field  processing.  We  will 
attempt  to  extend  matched  field  processing  to  larger  apertures:  a  200  km  aperture  (the  KNET)  and.  if  data  permit,  to  an  aperture  comprised  of  .several  Kazakh  arrays.  We  also  will 
investigate  the  potential  of  developing  matched  field  proce.ssing  to  roughly  locate  and  cla.ssify  natural  seismicity,  w  hich  is  more  diffuse  than  the  concentrated  sources  of  mining 
explosions  that  we  have  investigated  to  date. 
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ABSTRACT 

We  have  adapted  matched  field  processing — a  method  developed  in  underwater  acoustics  to  detect  and  locate 
targets—  to  classify  transient  seismic  signals  arising  from  mining  explosions.  Matched  field  processing,  as  we  apply 
it,  is  an  empirical  technique,  using  observations  of  historic  events  to  calibrate  the  amplitude  and  phase  structure  of 
wavefields  incident  upon  an  array  aperture  for  particular  repeating  sources.  The  objective  of  this  project  is  to 
determine  how  broadly  applicable  the  method  is  and  to  understand  the  phenomena  that  control  its  performance. 

We  obtained  our  original  results  in  distinguishing  events  from  ten  mines  in  the  Khibiny  and  Olenegorsk  mining 
districts  of  the  Kola  Peninsula,  for  which  we  had  exceptional  ground  truth  information.  In  a  cross-validation  test, 
some  98.2%  of  549  explosions  were  correctly  classified  by  originating  mine  using  just  the  Pn  observations  (2.5-12.5 
Hz)  on  the  ARCES  array  at  ranges  from  350  -  410  kilometers.  These  results  were  achieved  despite  the  fact  that  the 
mines  are  as  closely  spaced  as  3  kilometers.  Such  classification  performance  is  significantly  better  than  predicted  by 
the  Rayleigh  limit.  Scattering  phenomena  account  for  the  increased  resolution,  as  we  make  clear  in  an  analysis  of  the 
information  carrying  capacity  of  Pn  under  two  alternative  propagation  scenarios:  free-space  propagation  and 
propagation  with  realistic  (actually  measured)  spatial  covariance  structure.  The  increase  in  information  capacity  over 
a  wide  band  is  captured  by  the  matched  field  calibrations  and  used  to  separate  explosions  from  very  closely-spaced 
sources.  In  part,  the  improvement  occurs  because  the  calibrations  enable  coherent  processing  at  frequencies  above 
those  normally  considered  coherent. 

We  are  investigating  whether  similar  results  can  be  expected  in  different  regions,  with  apertures  of  increasing  scale 
and  for  diffuse  seismicity.  We  verified  similar  performance  with  the  closely-spaced  Zapolyami  mines,  though 
discovered  that  it  may  be  necessary  to  divide  event  populations  from  a  single  mine  into  identifiable  subpopulations. 
For  this  purpose,  we  perform  cluster  analysis  using  matched  field  statistics  calculated  on  pairs  of  individual  events  as 
a  distance  metric.  In  our  initial  work,  calibrations  were  derived  from  ensembles  of  events  ranging  in  number  to  more 
than  100.  We  are  considering  the  performance  now  of  matched  field  calibrations  derived  with  many  fewer  events 
(even,  as  mentioned,  individual  events).  Since  these  arc  high-variance  estimates,  we  are  testing  the  use  of 
cross-channel,  multitaper,  spectral  estimation  methods  to  reduce  the  variance  of  calibrations  and  detection  statistics 
derived  from  single-event  observations. 

To  test  the  applicability  of  the  technique  in  a  different  tectonic  region,  we  have  obtained  four  years  of  continuous  data 
from  4  Kazakh  arrays  and  are  extracting  large  numbers  of  event  segments.  Our  initial  results  using  132  mining 
explosions  recorded  by  the  Makanchi  array  are  similar  to  those  obtained  in  the  European  Arctic.  Matched  field 
processing  clearly  separates  the  explosions  from  three  closely-spaced  mines  located  approximately  400  kilometers 
from  the  array,  again  using  waveforms  in  a  band  (6-10  Hz)  normally  considered  incoherent  for  this  array.  Having 
reproduced  ARCES-type  performance  with  another  small  aperture  array,  we  have  two  additional  objectives  for 
matched  field  processing.  We  will  attempt  to  extend  matched  field  processing  to  larger  apertures:  a  200  km  aperture 
(the  KNET)  and,  if  data  permit,  to  an  aperture  comprised  of  several  Kazakh  arrays.  We  also  will  investigate  the 
potential  of  developing  matched  field  processing  to  roughly  locate  and  classify  natural  seismicity,  which  is  more 
diffuse  than  the  concentrated  sources  of  mining  explosions  that  we  have  investigated  to  date. 
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OBJECTIVE 

The  objectives  of  this  project  are 

•  To  investigate  the  limits  of  empirical  matched  field  processing  and  other  coherent  array  detection  and  parameter 
estimation  methods  as  receiver  aperture  size  increases  from  a  few  kilometers  to  many  hundreds  of  kilometers. 

•  To  investigate  techniques  for  extending  the  geographical  source-region  footprint  over  which  empirical  matched 
field  processing  and  other  coherent  calibrated  methods  apply. 

We  will  begin  by  reanalyzing  data  from  the  European  Arctic  in  order  to  reconfirm  the  potential  of  empirical  matched 
field  processing  that  has  been  indicated  by  our  work  under  a  previous  contract.  We  then  will  proceed  to  study  the 
central  Asia  region  to  assure  programmatic  relevance  and  to  exploit  the  large  belts  of  natural  and  man-made 
seismicity  to  test  the  applicability  of  the  technique  to  diffuse  seismicity. 

RESEARCH  ACCOMPLISHED 

The  effort  was  focused  initially  on  examining  the  degree  to  which  empirical  matched  field  processing  (EMPF)  can 
classify  signals  recorded  on  the  ARCES  array  in  northern  Norway  from  closely  spaced  mining  explosions.  There 
were  many  good  reasons  for  testing  the  new  algorithms  on  this  data  set.  Through  close  collaboration  with  scientists  at 
the  Kola  Regional  Seismological  Center  in  Apatity,  and  with  the  operators  of  numerous  mines  on  the  Kola  Peninsula, 
a  superb  database  of  Ground  Truth  information  has  been  accumulated  (Harris  et  al.,  2003),  consisting  of  many 
hundreds  of  events  from  a  period  of  several  years.  In  addition,  the  parameter  estimation  capabilities  of  the  ARCES 
array  for  such  signals  have  been  explored  in  depth.  The  ability  of  the  array  to  resolve  source  location  using  the 
classical  array-processing  paradigm  is  relatively  well  understood,  providing  an  excellent  basis  for  comparison  with 
the  new  technology.  The  panel  to  the  left  of  Figure  1  shows  the  location  of  the  sources  with  respect  to  ARCES. 


Figure  1.  Locations  of  sites  of  industrial  and  repeating  seismicity  in  relation  to  the  ARCES  array  (left)  and 
Google  Earth  image  of  the  Zapoljarni  open  cast  mining  region  with  the  ZPl  (Zapadny)  and  ZP2 
(Central)  mines  labelled.  The  mines  are  separated  by  approximately  2.5  km  and  are  approximately 
200  km  from  ARCES. 

The  new  data  from  Kazakhstan  and  Kyrgyzstan  has  allowed  the  application  of  the  same  techniques  to  a  different 
tectonic  setting,  of  programmatic  relevance,  with  high  levels  of  natural  and  artificial  seismicity.  Three  different 
apertures  of  sensor  configuration,  with  which  the  ability  to  expand  Coherent  Array  Processing  to  greater  sensor 
separations  can  be  examined,  are  present.  At  the  smallest  scale,  the  four  single  arrays  in  Kazakhstan  can  be  processed 
in  the  same  way  as  ARCES.  These  arrays  arc  however  somewhat  larger  and  far  more  sparse  than  ARCES  which  will 
have  consequences  for  waveform  semblance,  in  particular  at  the  higher  frequencies.  At  an  intermediate  scale,  the 
stations  of  KNET  constitute  a  network  with  an  aperture  of  the  order  200  km  and  a  relatively  broad  spectrum  of  sensor 
separations.  The  largest  aperture  to  be  considered  is  that  formed  by  the  four  Kazakh  arrays  as  a  single  network. 

Ground  Truth  of  the  kind  which  we  have  exploited  in  our  examination  of  industrial  seismicity  on  the  Kola  Peninsula 
is  not  yet  available  for  Kazakhstan.  This  will  complicate  the  evaluation  of  the  event  classification  attempted  under  the 
current  contract.  We  have  however,  been  assisted  greatly  by  access  to  classification  results  obtained  by  MacCarthy  et 
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al.  (2008)  and  Hartse  et  al.  (2008).  These  studies  are  aimed  at  event  sereening  and  souree  identifieation  using,  in 
partieular,  waveform  similarity  teehniques  and  satellite  imagery.  Their  original  souree  of  information  is  the  bulletin  of 
the  Kazakh  National  Data  Center  (KNDC). 

Deteeted  seismie  phases  are  usually  eharaeterized  by  the  parameters  Vapp  (apparent  veloeity)  and  0  (baekazimuth) 
for  the  purpose  of  phase  identifieation  and  assoeiation,  and  subsequent  event  loeation.  In  the  estimation  proeess,  these 
parameters  are  eoded  in  the  slowness  veetor,  s,  defined  by 

s  =  {sj:,Sy)  where  s*=ssin(©),  Sy  =  scos(0)  and  s  =  l/v^pp  (•) 


If  Xj  denotes  the  eoordinates  of  sensor y,  relative  to  the  array  referenee  site,  then  we  denote  the  time-series  reeorded 
at  this  site  by  r(t,  Xj). 

T 

For  a  given  frequeney  w  ,  the  steering  veetor  ff(a;,  s)  =  .  g-twS.XAf  j  ^Hows  a  measurement  of 

the  energy  ineident  on  the  array  eonsistent  with  the  plane  wavefront  hypothesis,  s,  using  (e.g.  Capon,  1969) 

2 


P{u;,3)  = 


■'  j=i 


=  ff(a;,  s)^ R{uj)€{(jj^  s) 


(2) 


where  denotes  the  number  of  sensors  in  the  array,  //denotes  the  Hermitian  transpose,  and  the  elements  of  the 
spatial  eovarianee  matrix,  H(a;),  are  given  by 


fimn(w)  =  r{t,Xm)e  r{t,Xn)e"^^dt^. 


(3) 


When  estimating  the  slowness  veetor,  s,  whieh  best  fits  the  incoming  wavefront,  it  is  typical  to  evaluate  a  normalized 
relative  beam-power  funetion  of  the  form 


5-  . 

= - E,<r(flK)> - 


(4) 


over  a  grid  of  slowness  space  and  summed  (incoherently)  over  a  wide  range  of  frequencies  j. 
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Figure  2.  Slowness  estimates  at  ARCES  using  broadband  f-k  analysis  for  first  regional  P-phases  from  the  ZPl 
(dark  blue)  and  ZP2  (light  blue)  mines.  The  true  baekazimuth  is  approximately  91.8  degrees  for  both 
mines  as  indicated  by  the  vertical  dashed  lines.  The  top  left  panel  displays  estimates  taken  from  the 
fully  automatic  ARCES  detection  lists  and  the  remaining  panels  arc  reprocessed  in  the  fixed 
frequency  bands  as  indicated.  The  41  events  from  the  ZPl  mine  and  the  25  events  from  the  ZP2  mine 
all  took  place  between  October  2001  and  September  2002. 
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The  right  panel  of  Figure  1  shows  a  satellite  image  of  two  mines,  denoted  ZPl  and  ZP2,  which  are  separated  by  less 
than  3  km.  Figure  2  shows  four  panels  of  slowness  estimates  for  initial  P-arrivals  from  these  two  mines  over  a  12 
month  period  from  Oetober  2001  to  September  2002. 

The  top  left  panel  (in  whieh  measurements  are  taken  from  the  fully  automatie  deteetion  lists)  shows  a  very  high 
spread  in  the  measured  baekazimuth  and  no  separation  between  the  two  populations.  (Under  standard  array 
proeessing  eonsiderations,  we  would  not  antieipate  being  able  to  resolve  these  two  nearby  sourees  with  a  single 
seismie  array  at  this  distanee.)  The  frequeney  band  seleeted  in  the  automatie  deteetion  lists  is  determined  by  an 
automatie  algorithm  on  a  ease-by-ease  basis  to  try  to  optimize  the  signal-to-noise  ratio  (SNR).  The  high  variability  of 
the  slowness  estimates  has  signifieant  eonsequenees  for  the  aeeuraey  of  fully-automatie  event  bulletins  (e.g.,  Ringdal 
and  Kvsema,  1989)  and  all  undertakings  whieh  ean  result  in  more  stable  parameter  estimates  should  be  investigated. 

The  three  remaining  panels  are  interesting  as  they  demonstrate  (a)  that  the  stability  is  very  mueh  improved  by 
seleeting  a  fixed  frequeney  band,  (b)  that  the  bias  and  spread  of  the  slowness  estimates  are  very  different  for  the 
different  frequeney  bands,  and  (e)  that  there  is  a  signifieant  separation  in  parameter  spaee  for  the  two  different  mines. 
Whilst  this  separation  is  evident  visually,  it  would  be  very  diffieult  to  identify  the  souree  eonfidently  just  given  the 
parameter  estimates  for  the  different  frequeney  bands.  It  is  the  intention  that,  in  bypassing  the  plane-wavefront 
parametrization,  that  Empirical  Matched  Field  Proeessing  will  be  able  to  exploit  the  eharaeteristie  souree  and  path 
effeets  responsible  for  the  patterns  shown  in  order  to  identify  a  likely  souree. 

If  we  calculate  spatial  eovarianee  matrices,  il(a;) ,  for  seismie  arrivals  known  to  have  originated  from  a  given  souree, 
we  ean  form  an  ensemble  eovarianee  matrix  from  a  speeified  population  a  of  events.  If  we  replaee  the 

theoretieal  plane  wavefront  steering  veetors  c(a;,  s)  in  Equation  4  with  the  prineipal  eigenveetors  of  Rq{^)9  we  ean 
obtain  an  equivalent  normalized  measure  of  the  energy  ineident  on  the  array  for  the  “a”  hypothesis: 

EitrlRiuj)} 


Ensemble  spatial  eovarianee  matriees  -Ra(^)  were  ealeulated  for  the  two  populations  of  events  displayed  in  Figure  2 
and  the  “matched  field  statistic”  specified  in  Equation  5  was  evaluated  for  both  these  empirieal  steering  veetors  for 
eaeh  event  known  to  have  taken  plaee  from  the  same  two  mines  during  the  period  2003-2004  (sec  Figure  3). 


Whilst  the  separation  of  the  ZPl  and  ZP2  populations  against 
the  ealibrated  empirieal  steering  vectors  is  reasonable,  many 
symbols  lie  in  a  zone  of  high  uneertainty  (shaded  grey)  where 
the  mateh  of  the  measured  wavefield  is  quite  similar  for  both 
of  the  souree  hypotheses.  Two  events  from  the  ZPl  mine 
have  been  miselassified  as  belonging  to  the  ZP2  eluster  (with 
a  marginal  differenee  in  deteetion  statisties),  although  we  do 
not  rule  out  errors  in  the  Ground  Truth  eolleetion. 

An  ideal  elassifieation  would  see  all  of  the  dark  blue  symbols 
in  the  lower  right  eomer  and  all  the  pale  blue  symbols  in  the 
upper  left  eomer. 


Figure  3.  Matched  field  statistics  as  defined  in  Equation  5  for  empirical  steering  vectors  calibrated  from  the 
2001-2002  populations  for  the  ZPl  and  ZP2  mines  on  the  x  and  y  axes  respectively.  Dark  and  light 
blue  symbols  indicate  that  the  spatial  covariance  matrices  were  calculated  from  signals  from  ZPl 
and  ZP2  events  respectively  during  the  period  2003-2004.  Frequeney  band  used:  2.5  -  12.5  Hz. 
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There  is  a  clear  need  to  understand  why  the  classification  performance  is  as  shown  in  Figure  3  and  investigate  which, 
if  any,  improvements  are  possible.  The  ability  to  attribute  candidate  signals  to  one  or  other  of  the  sources  will  depend 
upon  how  well  the  principal  eigenvector  of  the  ensemble  spatial  covariance  matrix  for  each  cluster  represents  the 
spatial  covariance  observed  for  individual  events  associated  with  that  source,  and  the  degree  to  which  the  empirical 
steering  vectors  from  the  different  sources  differ  from  each  other.  To  this  end,  we  perform  cluster  analysis  upon 
empirical  steering  vectors  calculated  from  individual  events.  Specifically,  we  calculate  a  similarity  measurement  as 
defined  in  Equation  6  between  the  principal  eigenvectors  of  the  spatial  covariance  matrices  between  any  two  distinct 
data  segments.  Note  that  the  measure  is  an  average  value  taken  over  a  wide  range  of  frequency  bands.  Figures  4  and  5 
show  cross-matrices  and  corresponding  single-linkage  dendrograms  for  the  ZPl  and  ZP2  event  sets  for  2001-2002. 
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Figure  4.  Similarity  matrix  and  corresponding  single-linkage  dendrogram  for  events  from  the  ZPl  mine 
2001-2002. 
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Figure  5.  Similarity  matrix  and  corresponding  single-linkage  dendrogram  for  events  from  the  ZP2  mine 
2001-2002. 

It  is  clear  from  Figure  5  that  the  ensemble  of  events  selected  to  represent  the  ZP2  mine  consists  of  two  distinct 
clusters  according  to  the  similarity  metric  defined  in  Equation  6.  We  then  repeat  the  classification  exercise  displayed 
in  Figure  3  but  replacing  the  ZP2  matched  field  statistic  with  the  maximum  of  two  different  statistics,  calculated  from 
the  ensemble  covariance  matrices  from  each  of  the  two  clusters  evident  in  Figure  5.  It  is  clear  in  the  right  hand  panel 
of  Figure  6  that  the  classification  is  far  better  than  that  displayed  in  Figure  3  where  only  a  single  empirical  steering 
vector  was  available  to  represent  the  whole  ZP2  source  region.  Almost  no  events  appear  in  the  grey  shaded  region 
when  the  two  steering  vector  representation  was  used  for  the  ZP2  site.  The  same  two  events  in  the  ZPl  population 
remain  miselassified,  and  there  is  evidence  in  Figure  4  for  (at  least)  two  events  which  appear  to  match  quite  poorly. 
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Figure  6.  Waveforms  from  the  ZP2  events  (left)  with  traces  colored  according  to  the  classification  displayed  in 
Figure  5,  and  Matched  Field  Statistics  (right)  with  ZP2  events  represented  by  two  different  empirical 
steering  vectors.  The  empirical  steering  vectors  for  both  mine  clusters  are  calculated  from  the 
2001-2002  event  pools  and  the  dark  and  light  symbols  represent  ZPl  and  ZP2  events  respectively 
from  2003-2004.  Frequency  band  used:  2.5  -  12.5  Hz. 

We  conclude  that  cluster  analysis  is  necessary  as  a  prerequisite  to  the  calculation  of  ensemble  covariance  matrices  in 
order  that  the  validity  of  the  selected  cluster  can  be  determined.  If  a  cluster  clearly  separates  into  two  populations,  it  is 
likely  that  matched  field  statistics  calculated  using  an  ensemble  spatial  covariance  matrix  generated  from  the  full 
population  are  not  going  to  provide  a  good  classification.  Events  which  are  spurious  and  which  clearly  do  not  belong 
to  the  cluster  will  be  identified  using  the  procedure  described  here. 


The  scattering  effects  which  are  so  problematic  for  coherent  array  processing  of  high  frequency  phases  on  ARCES 
are  even  more  so  for  the  arrays  in  Kazakhstan  due  to  the  sparse  geometries  (see  Figure  7).  Whilst  the  problem  of  low 
SNR  at  low  frequencies  is  less  acute  than  at  ARCES,  the  SNR  still  improves  greatly  at  higher  frequencies  (sec 
Figure  8).  Estimating  the  slowness  using  broadband  f-k  analysis  is  a  trade-off  between  SNR  and  signal  coherence. 
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Figure  7.  Comparison  of  the  eonfigurations  of  the  MKAR  and  ARCES  arrays.  Whilst  both  arrays  have  the 

same  fundamental  design  (of  eoneentrie  rings  with  inereasing  numbers  of  sensors),  the  sparsity  of  the 
MKAR  array  makes  proeessing  of  regional  phases  very  diffieult.  The  inner-most  sensors  of  ARCES 
whieh  are  so  erueial  to  the  detection  and  estimation  of  high  frequeney  phases  are  missing  from 
MKAR. 
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The  data  segment  in  Figure  8  is  the  result  of  a 
mining  explosion  on  June  30,  2007,  at  a 
distance  of  approximately  400  km  from  the 
MKAR  array.  This  event  belongs  to  the  cluster 
labelled  344  by  Hartse  et  al.  (2008). 

The  SNR  for  the  Pn  arrival  is  high,  but  only 
for  frequencies  above  4  Hz.  This  is  unfortunate 
from  the  point  of  view  of  phase  estimation 
using  traditional  array  processing  since  (see 
Figure  9)  almost  all  slowness  estimates  above 
4  Hz  are  qualitatively  incorrect  due  to 
scattering  and  incoherence  together  with 
sidelobes  in  the  array  response  function  for 
these  frequencies.  Many  of  the  estimates 
below  4  Hz  arc  qualitatively  misleading  due  to 
a  low  SNR. 
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Figure  8.  A  three-minute  long  data  segment  on  the  single  channel  MK01_SHZ  filtered  in  various  frequency 
bands  as  indicated.  Sec  above  text  for  details. 
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Figure  9.  The  left  hand  panels  show  scatter-plots  of  slowness  estimates  in  the  frequency  bands  indicated  for 
initial  Pn  arrivals  at  MKAR  for  events  in  clusters  310,  328,  and  344  in  Hartse  et  al.  (2008).  The 
anticipated  backazimuth  for  all  events  is  close  to  340^  and  the  apparent  velocity  should  be  close  to  8 
kms  ^  The  f-k  grids  to  the  right  display  the  relative  beam-power  as  a  function  of  slow  ness  for  the 
arrival  at  MKAR  at  a  time  2007-181:06.30.31.425  (see  Figure  8). 
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Figure  9  indicates  at  a  glance  that  an  attempt  at  source  classification  based  upon  direction  estimates  from  traditional 
plane  wavefront  array  processing  are  not  likely  to  succeed  given  the  abysmal  coherence  on  the  MKAR  array  for  high 
frequency  regional  phases.  We  attempt  a  classification  using  empirical  matched  field  processing;  given  that  the 
clusters  were  derived  by  Hartse  et  al.  (2008)  and  MacCarthy  et  al,  (2007)  using  waveform  correlation  methods,  some 
degree  of  success  is  anticipated.  Cluster  analysis  was  performed  on  the  empirical  steering  vectors  for  each  of  the 
individual  events  prior  to  an  attempt  at  classification  using  ensemble  covariance  matrices.  The  results  are  displayed  in 
Figure  10.  From  this  plot,  it  is  clear  that  multiple  empirical  steering  vectors  are  necessary  to  characterize  the 
wavefields  recorded  by  the  range  of  events  from  each  of  these  clusters. 


Figure  10.  Matched  field  similarity  matrix  compiled  over  three  distinct  clusters  (310,  328,  and  344)  from  the 
study  of  Hartse  et  al.  (2008).  A  number  of  distinct  subpopulations  arc  revealed  from  the  eluster 
analysis  of  the  empirical  steering  vectors  for  the  Pn  arrival  at  the  MKAR  array  (espeeially  for 
cluster  328)  and  these  need  to  be  taken  into  eonsideration  prior  to  matehed  field  elassifieation  using 
ensemble  spatial  covariance  matrices. 

Figure  1 1  shows  that  the  matched  field  classification  has  attributed  each  of  the  individual  events  to  the  correct  source. 
We  have  therefore  demonstrated  that  the  spatial  covariance  matrices  provide  a  coherent  array-processing  technique 
which  works  at  high  frequencies  over  an  aperture  for  which  simply  delay-and-stack  beamforming  and  broadband  f-k 
analysis  fail. 
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Figure  1 1 .  Classiflcation  results  for  clusters  3 1 0, 328,  and  344,  using  multiple  ensemble  covariance  matrices  for 
each  of  the  subpopulations  displayed  in  Figure  11. 

CONCLUSIONS  AND  RECOMMENDATIONS 

We  have  demonstrated  that  empirieal  matched  field  processing  is  able  to  elassify,  with  a  high  degree  of  confidence, 
events  from  sources  of  repeating  seismicity  even  for  cases  where  classical  waveform  correlation  methods  become 
problematic  due  to  differences  in  the  source  time  function.  This  has  also  been  demonstrated  for  an  array  with  an 
aperture  over  which  incoherence  defeats  traditional  methods  of  array  processing. 

One  of  the  most  important  aspects  of  the  procedure  that  has  come  to  light  during  this  reporting  period  is  the  necessity 
of  performing  cluster  analysis  on  the  empirical  steering  vectors  prior  to  the  building  of  ensemble  covariance  matrices. 
If  an  attempt  is  made  to  use  an  ensemble  covariance  matrix  constructed  from  significantly  differing  subpopulations, 
this  will  result  in  empirical  steering  vectors  which  represent  the  entire  population  relatively  poorly. 

It  is  of  great  interest  as  to  how  many  events  arc  required  to  characterize  a  specified  source  region  adequately;  there 
arc  clearly  circumstances  where  one  would  hope  to  achieve  a  reasonable  performance  with  only  a  single  observation. 
Part  of  the  motivation  for  the  construction  of  the  ensemble  covariance  matrices  is  to  reduce  the  variance  in  the 
narrow-band  spatial  covariance  matrices  (Equation  3).  Multitaper  methods  (Thomson,  1982)  produce  spectral 
estimates  with  low  variance  and  the  recently  published  codes  of  Prieto  ct  al.  (2009)  include  routines  for  robust  and 
stable  multitaper  coherence  estimation.  We  have  implemented  these  routines  into  our  codes  and  are  currently 
investigating  the  variability  of  the  estimates  of  the  spatial  covariance  matrices. 
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